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ABSTRACT 

Visual and ultraviolet spatially resolved (~ O'.'l) spectra of SN 1987A 
obtained on days 3715 and 3743 with the Space Telescope Imaging Spectrograph 
on the Hubble Space Telescope show that the high- velocity SN debris is colliding 
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with circumstellar gas. Very broad Lya emission with velocities extending to 
~ ±20, 000 km s _1 originates inside the inner circumstellar ring and appears 
to fill most of the surface area within 0V67 ± / .'03 (0.14 pc at a distance of 
50 kpc) of the ring's center. The observed Lya flux from the shocked ejecta 
is (1.85 ±0.53) x 10~ 13 erg cm" 2 s' 1 and (1.25 ± 0.51) x 10~ 12 erg cm" 2 s" 1 
after correcting for extinction. A spatially unresolved blue-shifted emission 
feature was discovered in Ha (and other lines) on the inner ring at p. a. 31° ± 8°. 
The Ha emission extends to —250 km s -1 with no corresponding red-shifted 
emission. This highly localized interaction appears to be the initial contact of 
the supernova blast wave with an inward protrusion of the inner ring. The broad 
Lya emission and the 'hot spot' are separate interaction phemonena associated 
with the reverse and forward shocks, respectively. We also find that the size 
of the inner ring in forbidden lines of oxygen has a dependence on ionization 
potential, in agreement with photoionization models of the ring. 

Subject headings: supernovae: individual(SN 1987A) - supernova remnants - 
circumstellar matter 



1. INTRODUCTION 

SN 1987A in the Large Magellanic Cloud continues to provide new opportunities to 
study the late-time evolution of supernovae and to probe their stellar origins. The collision 
of the ejecta with the progenitor's circumstellar (CS) medium has been anticipated since 
its discovery in mid-1987, due to its proximity to the supernova (Fransson et al. 1989). The 
coincident emergence of radio (Staveley-Smith et al. 1992, 1993) and X-ray (Beuermann 
et al. 1994) emission in 1990, their subsequent monotonic increase, and the deceleration 
of the supernova shock (Staveley-Smith et al. 1993) indicate that the high-velocity SN 
shock encountered low-density CS gas starting around 1200 days after the explosion. This 
gas corresponds to mass lost by the progenitor near the end of its red supergiant (RSG) 
phase (Chevalier & Dwarkadas 1995, hereafter CD95). Gaensler et al. (1997) find the radio 
emitting region is located just inside the central ring, based on a ~0'/4 resolution 9 GHz 
map, and is expanding at ~ 2,800 km s _1 (12 mas yr _1 ). The shock interaction region, 
bounded by forward and reverse shocks, consists of an outer zone of shocked CS gas and an 
inner zone of shocked ejecta (Chevalier 1982). 

The SN 1987A CS medium is a fossil record of the pre-explosion mass loss and may 
hold important clues to the nature of the progenitor (e.g. single vs. binary star). The gas in 
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the central ring (n e ~ 10 4 cm~ 3 ) is significantly enhanced in helium and nitrogen, indicating 
that the progenitor was in a post He-core burning phase when it exploded (Fransson et al. 
1989, Sonneborn et al. 1997). The presence of lower density gas (n e ~ 100 cm -3 , CD95) 
inside the ring traces the mass loss at the end of the RSG stage, after formation of the 
inner ring. 

The small angular extent and complex structure of the SN 1987A CS gas and debris 
require the high spatial resolution of the Hubble Space Telescope (HST). WFPC2 and 
FOC imagery revealed the remarkable thin-ring structures (Jakobsen et al. 1991, [O III]; 
Burrows et al. 1995, Ha; Plait et al. 1995, [N II]). However, the first generation HST 
spectrographs had one-dimensional detectors and could only provide spatial information 
on the scale of the available apertures. The Space Telescope Imaging Spectrograph (STIS) 
provides imaging spectroscopy by combining excellent optical performance of imaging 
optics and two-dimensional detectors (Kimble et al. 1997). With good spectral resolution 
and i/ST-limited spatial resolution, we are now able to study the spatial dependence of 
the line emission from the inner ring and do the first spectral analysis of the interaction of 
the SN ejecta with the CS gas. In an accompanying paper, Michael et al. (1997, hereafter 
M97) model the CS interaction and the emission processes that produce the Lja emission 
discussed below. 

2. OBSERVATIONS 

2.1. Ultraviolet Spectral Imaging 

Ultraviolet spectral images of SN 1987A and the inner CS ring were taken with the 
STIS FUV-MAMA detector on 1997 May 24 (3743 days after core collapse). Two exposures 
(Pgm ID 7123) were taken in time-tag mode with the G140L grating (~ 1130 - 1720 A) for 
a total exposure time of 4180 sec. The plate scale of the image is 0'.'0244 per pixel in both 
the spatial and spectral directions. The spatial resolution is ~ 0'.'06. The G140L spectral 
resolution is ~ 1.2 A (0.5833 A/pixel). The 2" x 2" aperture was chosen to include the 
entire inner circumstellar ring (1'.'7 x 1'.'2, Plait et al. 1995) and exclude the neighboring B 
stars, primarily star 3 (V ~ 16, l'/6 from the SN; see Fig. la). The G140L spectral image is 
shown in Figure 2. 

A striking and previously unobserved emission feature, extending from <; 1145 to 
~ 1300 A, is identified as Lya from the shock interaction of the high- velocity debris with 
CS gas. The strong radiation filling the square aperture is geocoronal Ly<x The diffuse 
continuum extending the length of the spectrum is UV background from the LMC. The 
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detector background count rate is extremely low (< 1 x 10 5 cts s 1 px 1 , Kimble et 
al. 1997) and is negligible in this data. Images of the inner ring are also detected in 
O IV] AA1397-1407, N IV] A1486.5, and He II A1640.5; their orientation indicates that the 
aperture's spatial axis is at p. a. 115° ±4°. A weak, spatially extended UV continuum (1400 
- 1700 A) from the debris is also present. 

Wavelength calibration exposures using the Pt-Cr/Ne line lamp were taken through 
the 52" x O'.'l slit for each G140L exposure. The absolute wavelength scale was referenced 
to the small offset (Ax) of the SN from the center of the aperture. We find the SN offset to 
be Ax = +0"12 ± 0'.'15 in the dispersion direction and Ay = — 0'!17 ± 0'.'04 in the spatial 
direction, as determined from the center of the ring in O IV] , N IV] , and He II . 

2.2. Visible Spectral Imaging 

Medium-resolution spectral images (Pgm ID 7123) of the inner ring were obtained in 
two grating settings on 1997 April 26 (day 3715): G750M (6295 - 6867 A, resolution = 
1.02 A, 643 sec), and G430M (4818 - 5104 A, resolution = 0.83 A, 828 sec). The 52" x 2" slit 
was chosen to include the entire ring and was oriented along its major axis (Fig. lb), thereby 
providing maximum separation in the dispersion direction between images of the ring in 
adjacent emission lines. The narrow emission lines from the ring (vfwhm ^ 15 km s _1 ) 
remain spectrally unresolved in these observing modes (Av ~ 50 km s _1 ), while the debris 
spectral features (wfwhm ~ 2500 km s^ 1 ) are highly dispersed. In addition to the CS ring 
and the SN debris, the G750M and G430M spectral images (see Figure 3) also show diffuse 
Ha, H/3, and [O III] emission filling the long slit and the continuous spectra of star 3 
(including stellar Ha and H./3 emission lines) and several fainter stars. 

An [O II] A3727 image (Pgm ID 7122) of SN 1987A and its ring system was taken 
with STIS on 1997 April 4 (day 3693), using the F28X50OII filter. Four exposures totaling 
2415 sec were obtained. 

3. RESULTS 

3.1. Ultraviolet Spectroscopy 

The Lya emission from SN 1987A reported here is the signature of the shock interaction 
between high-velocity SN ejecta and CS gas inside the bright inner ring. Two aspects of the 
Lya feature are particularly noteworthy: a) the spatial and velocity distribution of the Lya 
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flux, and b) the relative brightnesses of the red and blue wings of the line profile. 

We measure the red edge of the Lya profile at 1301.7 ± 3.7 A, or +20, 900 ± 900 km s" 1 
relative to 1216.825 A, Lya in the LMC rest frame. The red edge could, however, be 
as large as 1328 A (27,400 km s" 1 ), depending on the details of the diffuse background 
subtraction. A Doppler shift of —20, 900 km s" 1 on the blue side of the profile occurs 
at 1131 A. Unfortunately, we cannot measure the profile shortward of ~ 1145 A due to 
decreasing ifST/STIS sensitivity below 1170 A. Measurement of the blue wing of the Lya 
profile awaits observations by the Far Ultraviolet Spectroscopic Explorer (FUSE) in 1999. 

The Lya brightness profile in the spatial direction (across the spectrum) sets a limit 
on the spatial extent of the shock interaction region in a direction close to the major axis of 
the ring. Figure 4 shows spatial brightness profiles in the blue and red wings of Lya. The 
Lya surface brightness is fairly uniform and there is no significant brightening along the 
upper and lower edges of the feature (cf. M97). The spatial extent is 1'.'22 ±0'.'06 and 1'.'14 
± O'.'ll on the red (1250 ± 8 A) and blue (1179 ± 8 A) sides of the aperture. The projected 
size of the ring in the observing geometry is 1'.'51 ± 0"06, based on the ring size in [O III] 
(Plait et al. 1995). Thus the interaction region is located at 79% ± 4% of the ring's radius, 
or ~ 5 x 10 17 cm from the center of the debris. The spatial extent of the interaction region 
in the dispersion direction cannot be measured directly from the STIS data, due to the 
fact that the observed (x, y) distribution of Lya photons is a convolution of the geometry 
and the velocity field. The separation of the two effects is model-dependent (M97), but 
the spatial extent is probably ^ 1'.'4. The blue and red wings are tapered toward higher 
velocities and the maximum Doppler shift is roughly aligned with the center of the ring. 

The Lya line profile (Fig. 5) is the sum of the observed count rate at each wavelength 
point multiplied by the G140L sensitivity function. The integrated observed flux assumes 
that the part of the Lya profile obscured by geocoronal emission decreases linearly 
from 1187.0 A to 1241.3 A (see Fig. 5), corresponding to the edges of the 2"-wide 
aperture in the G140L image. The points where the profile meets the background are 
1144.4 and 1301.7 A. The integrated Lya flux is (1.85 ± 0.53) x 10~ 13 erg cm' 2 s _1 , 
consistent with the hydrodynamical models of Borkowski, Blondin, & McCray (1997). The 
dereddened Lya flux is (1.25 ± 0.51) x 10~ 12 erg cm" 2 s _1 (E(S — V) — 0.16, Fitzpatrick 
& Walborn 1990). The dereddened fluxes from the observable portions of the profile are 
(3.3 ±0.9) x 10~ 13 erg cm" 2 s" 1 (1144.4 - 1187.0 A) and (2.3 ± 0.5) x 10~ 13 erg cm" 2 s" 1 
(1241.3 - 1301.7 A). This asymmetry is probably due to the ~10-month light travel time 
across the interaction region, implying a short time scale for changes in the shock spectrum. 
Dust in the ejecta might also be a factor. There is no conclusive evidence for N V A1240 
emission from the shocked gas (see M97). The intensity variations in the red wing may be 
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the result of strong interstellar absorption lines (e.g., 1260, 1300, and 1335 A, cf. Sonneborn 
et al. 1997) convolved with the spatial extent of the interaction region. 

There are significant uncertainties in the Lya flux measurements. The net effect of 
statistical errors is very small (< 1%) and is negligible by comparison with the other error 
sources: the geometry of the interaction region, systematic errors in the flux measurement, 
interstellar extinction, and the HST /STIS G140L sensitivity function. 

We estimate that the uncertainty in the evaluation of the Lya flux hidden by 
geocoronal emission (see Fig. 5) is ~25%. The estimated uncertainty in the preliminary 
G140L sensitivity function used for this flux calculation is ~ 20% longward of ~ 1240 A 
and increases to ~ 40% at the shortest wavelengths. (The post-launch STIS calibration 
is currently underway and the expected uncertainty for the complete sensitivity function 
should be ^ 10% for all wavelengths). Other potential sources of systematic error include 
the selection of the end points of the integration regions and background subtraction. 

The spatial extent of the CS interaction region in the dispersion direction 
(±0'.'7 = ±16.7 A = ±4100 km s _1 ) introduces uncertainties in the calibrated Lya profile, 
unless its geometry is known and modelled. If the Lya surface brightness is symmetric, 
as may be the case for SN 1987A (M97), this effect would be significantly reduced. 
Nevertheless, we estimate this uncertainty to be 20 - 25% of the Lya flux. 

Errors in the Galactic and LMC extinction laws, in particular at A < 1200 A, are an 
additional source of error in the dereddened flux. The uncertainty in the extinction law near 
1200 A is about 10% (Fitzpatrick & Walborn 1990), while for 1140 - 1200 A we estimate 
it to be ~ 20%. The uncertainty in E(B-V) for SN 1987A is 0.02 mag (Fitzpatrick 
& Walborn 1990). The net uncertainty due to extinction is ~ 25%. Combining these error 
sources, the total estimated uncertainty in the observed Lya flux is 30 - 35% (38 - 43% for 
the dereddened flux). 

Apart from the Lya emission, the inner ring is detected in the G140L image 
in O IV] AA1397-1407, N IV] A1486.5, and He II A1640.5 (see Fig. 2), indicating 
that a significant fraction of hot gas remains in the ring. The dereddened surface 
brightnesses for O IV], N IV], and He II are (6.0 ± 1.4), (10.8 ± 2.3), and 
(15.4 ± 3.4) x 10~ 15 erg cm -2 s -1 arcsec -2 , respectively. In addition, a faint and 
spatially extended far-UV continuum of the SN debris (FWZI ~ 0'.'44 and centered on the 
O IV], N IV], and He II images of the ring) is detected longward of 1400 A with an average 
dereddened flux of (1.6 ± 0.6) x 10~ 16 erg cur 2 s _1 A -1 . Although a V ~ 20 A5 star (Plait 
et al. 1995) is superposed on the ring at p. a. 227° (see Fig. 1), such a star would contribute 
less than 5% of the observed flux. There is no evidence for a point source spectrum in the 
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G140L data. 



3.2. Visible Spectroscopy 

The STIS G430M and G750M observations of SN 1987A yield nearly monochromatic 
images of the inner CS ring in a large number of lines: H/3, [O III] AA4958.9, 5006.8, [O I] 
AA6300.3, 6363.8, [N II] AA6548.0, 6583.4, Ha, and [S II] AA6716.5, 6730.8, are observed. 
Parts of the outer rings are noticeable in Ha, [N II], and [O III]. Ha emission from the SN 
debris is detected as a very broad, spatially resolved emission feature passing through the 
center of the [N II] and Ha ring images. The observed images of the ring (Fig. 6) have 
similar morphology and intensity variations. 

These spectral images revealed a spatially unresolved, and previously unknown, blue- 
shifted Ha emission feature (also detected in H/3, [O I] A6300, and [O III] A5007) located 
close to or on the inner ring at p. a. 31° ± 8° (Pun et al. 1997). The Ha emission extends 
to —250 km s" 1 with no corresponding red-shifted emission (see Fig. 6d). Garnavich, 
Kirshner, & Challis (1997) pinpointed this 'hot spot' when they found a brightening at 
the same p. a., but located slightly inside the ring, in 1997 July WFPC2 Ha, [O III], and 
broad-band images. The hot spot must be a separate phenomenon from that producing 
the spatially-extended, high-velocity Lya emission. It may be a manifestation of the initial 
encounter of the SN blast wave with an inward protrusion of the inner ring. If so, we would 
expect it to continue to brighten rapidly and that more such spots will appear during the 
next few years. 

The STIS observations also help determine the location of the outer rings relative to 
the SN. A portion of the northern outer ring is projected on the SN debris (cf. Fig. 1). 
In Fig. 6d, the image of this ring in [N II] A6583 has a ~ 0'.'3 gap where it crosses the SN 
position, indicating that this section is located behind the debris. 

The size of the ring as a function of p.a. in [O I] A6300, [O II] A3727, [O III] A5007, and 
[S II] A6731 is shown in Figure 7. In most places the ring size in the oxygen lines increases 
with ionization stage. The radius is most reliably measured at the E ansa (p.a. 88°), 
where 9 = / . / 792 and 0'.'845 for [O I] and [O III], respectively. The [S II] A6716/A6731 and 
[N II] A6548/A6583 line ratios (Fig. 8) are more uniform than the individual line images 
(Fig. 6). The mean [S II] ratio is 0.54 ± 0.02, corresponding to N e = 6000 cm -3 . This 
is similar to that found by Lundqvist et al. (1997), based on an analysis of HST/FOS 
ring spectra from earlier epochs at p.a. ~300°. The presence of a significant ionization 
gradient in the ring and a fairly uniform density for the emitting gas are consistent with 
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the photoionization/recombination models (Lundqvist & Fransson 1996; Lundqvist & 
Sonneborn 1997). 

4. CONCLUSIONS 

STIS UV spectral imaging has further characterized the shock interaction region 
around SN 1987A. Broad Lyct emission originates inside the inner ring, at ~ 80% of the 
ring's radius, corresponding to the location of the reverse shock front (M97). The fairly 
uniform Lyct surface brightness and large Doppler shifts indicate that the shock geometry 
is not confined to the plane of the inner ring and may in fact be spheroidal, as shown 
by M97. This may provide evidence about the geometry of the RSG wind above and 
below the ring plane. The emerging picture of the debris-CS gas interaction region is close 
indeed to that proposed by CD95. The asymmetry between the blue and red wings of 
Lya indicates that the emission from the reverse shock may be evolving on a time scale of 
months. The location of the reverse shock is consistent with unimpeded debris expansion 
at ~ 40, 000 km s _1 up to ~ 1200 days, when the radio emission picked up, followed by 
expansion of the shock at ~ 2, 800 km s _1 to the present. 

The hot spot discovered at p. a. 31° appears to be a different, highly localized shock 
interaction phenomenon closer to the ring than that observed in Ly<x The location of the 
hot spot indicates that the forward shock has reached the innermost edge, or an inward 
protrusion, of dense ring material. This places the forward shock at ~ 1.1 times the the 
reverse shock radius. If this picture is correct, similar hot spots may light up around the 
ring in the coming years, providing a detailed map of the ring's radial geometry. 

We are grateful to Bob Kirshner and the SINS collaboration for sharing their recent 
WFPC2 results with us before publication and to Roger Chevalier and Claes Fransson for 
useful discussions. We thank the staff of the Space Telescope Science Institute and the 
STIS Team at GSFC and Ball Aerospace for their support and the HST Project for the 
opportunity to make these observations as part of the Early Release Observation program. 
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Figure Legends 

FIGURE I. - - The size and location of the STIS apertures are shown on a 1995 
composite WFPC2 image of SN 1987A (courtesy of P. Challis, R. Kirshner, and the SINS 
collaboration). The image size is 13" square, a) The 2"x 2" aperture location for the G140L 
observation, b) The 52" x 2" aperture orientation (p. a. 87°) for the G430M and G750M 
exposures. 

FIGURE 2. — a) STIS FUV-MAMA G140L (1130 - 1720 A) raw spectral image of 
SN 1987A. The image was reconstructed from time-tag mode data. The displayed image 
is a 1024 x 300 pixel subarray of the 1024 x 1024 image. Although star 3 lies outside 
the aperture, some light from it is still detected in the G140L image as the continuum 
emission along the top edge of the spectrum, b) Expanded view of Fig. 2a, with background 
corrections. Flux from star 3 was subtracted with row-by-row continuum fits to the data. 
The diffuse background was approximated by extracting the spectrum in the ~ 0"2 gap 
between the extended Lya emission and the light from star 3. This spectrum was subtracted 
from every row of the image after it was scaled to match the gradient in the background 
along the slit. Both images are displayed on a logarithmic intensity scale. Wavelength 
scales are shown in the LMC rest frame. 

FIGURE 3. — a) STIS G430M 485lA spectral image of SN 1987A shows the images of 
the inner circumstellar ring in H/9 and [O III] A4959 and A5007. These lines are also present 
in diffuse emission from the LMC filling the 52" x 2" aperture. The bright point source 
spectrum is from star 3. The image shown is a 950 x 450 pixel subarray of the 1024 x 1024 
image. The image was processed to remove cosmic rays and hot pixels, b) same as Fig. 3a, 
for the G750M 658lA exposure. The images are shown on a logarithmic intensity scale and 
the contrast is stretched to enhance weaker features. The spatial scale in the dispersion 
and cross-dispersion directions are not identical in the raw G430M and G750M images due 
to anamorphic distortion, a normal characteristic of gratings used in the off-Littrow mode. 
The corrected images shown here are stretched in the dispersion direction by factors of 
1.0776 (G750M/6581 A) and 1.1223 (G430M/4961 A) to achieve 0"0507 x 0"0507 pixels. 
Wavelength scales are shown in the LMC rest frame. 

FIGURE 4. — Spatial distribution of the Lya emission along the aperture. The spatial 
scale goes from the lower edge (0") to the top edge (2") . The cross-cut in the red wing was 
made at 1250 ± 8 A and at 1179 ± 8 A in the blue wing. The vertical dashed lines mark the 
location of the inner edge of the CS ring. 

FIGURE 5. — Calibrated STIS G140L spectrum of SN 1987A. The thick line indicates 
where the Lya profile was measured for the integrated flux. The CS Lya emission within 
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the aperture is approximated by the broken line. The feature at 1470 A is a combination of 
N IV] and SN debris. 

FIGURE 6. — STIS images of the SN 1987A inner CS ring, a) [S II] A6716 and A6731. 
b) [O II] A3727. c) [O III] A5007. d) [N II] A6548, Ha, [N II] A6583. e) H/3. f) [O I] A6300. 
Images a), c), d), e), and f) are subarrays of the images shown in Fig. 3. Wavelength 
increases from left to right. The [O II] direct image (b) was taken with the F28X50OII 
filter; the SN 1987A debris and star 3 are visible (cf. Fig. 1). The pixel size of these images 
is 0"05 x 0'.'05. 

FIGURE 7. — Size of the SN 1987A inner CS ring as a function of position angle in 
the light of several emission lines. 

FIGURE 8. - - Brightness ratios [S II] A6716/A6731 (right) and [N II] A6548/A6583 
(left) for the SN 1987A inner ring from the G750M exposure. For the [N II] ratio, only 
points with errors < 0.25 are included, a total of 373 points. 
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